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Abstract 

A Su-Schrieffer-Heeger model with added PT-symmetric boundary term is studied in 
the framework of pseudo-hermitian quantum mechanics. For two carefully chosen special 
cases, a complete set of pseudometrics is constructed in closed form. When complemented 
with a condition of positivity, the pseudometrics determine all the physical inner products 
of the considered model. 


1 Introduction 


The Su-Schrieffer-Heeger (SSH) model |1, 2] is one of the benchmark topologically nontrivial 
models in physics of condensed matter, ft provides a convenient description of certain physical 
systems [3, 4], and serves as one of the simplest examples of topological insulators [5]. The 
Ti-site SSH model may be expressed as 


N 

h ssh = [*4 - A cos e ) a li-i a 2 i + t(l + A cos 9)c 4a 2i+1 + h.c. 

2—1 


( 1 ) 


with a \, a, being the i-th site fermionic creation and annihilation operators. It has been 
suggested [2] to complement this hermitian SSH model with a non-hermitian TT-symmetric 
(invariant under simultaneous action of parity and time reversal) boundary term 


Hr = 7 a ! a i + l*°la n 
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( 2 ) 
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with 7 = ia, a G M. Despite apparent non-hermicity of the resulting operator H = 
Hssh + Hi, it has been shown that its eigenvalues remain real for |a| sufficiently small. This 
motivates further examination of this model from the viewpoint of 7T-,symmetric quantum 
mechanics. 

IPT-symmetric quantum mechanics [ 6 , 7, 8 ] is a theoretical framework for studying rep¬ 
resentations of quantum observables by non-hermitian operators. It is based on a a sim¬ 
ple fact, that an observable may be represented by a non-hermitian operator on a given 
Hilbert space 77 (F) = (V, (•]•)), as long as such operator is hermitian in another Hilbert space 
77 ■ s - ) = (V, {(fi\ 0)q) [9]. The inner products on V are in one-to-one correspondence with so-called 
metric operators 0 through the formula 

(0|'0)e = (01©!^) (3) 

where we require 0 to be bounded, hermitian and positive in order to generate a genuine 
inner product [10]. The hermicity condition of H in 77could be expressed in operator form as 

H^O = OH (4) 

which is a direct generalization of the standard condition H = 770 Operators H satisfying 
eq. 4 for at least one metric operator 0 are called quasi-hermitian [11], It is clear that on a 
given Hilbert space 77 (i y all admissible observables are described by quasi-hermitian instead of 
hermitian operators. Furthermore, recall that for any bounded, positive 0, there always exists 
a decomposition 0 = with H bounded. Inserting this decomposition into eq. 4 yields 

h := fiTTTT 1 = (H77fr 1 ) t = 70 (5) 

which shows that being quasi-hermitian is equivalent to being boundedly diagonalizable with 
real spectrum (instead of unitarily diagonalizable). However, the construction of the similarity 
transformation D may be in general very difficult for a given H. Special attention is usually 
given to parametric quasi-hermitian models with nontrivial domains of observability. The 
boundaries of such domains are formed of the so-called exceptional points (EPs) [12, 13], points 
in parameter space, where the operator ceases to be diagonalizable. The boundary crossings of 
physical parameters (e.g. time) are studied under the name quantum catastrophes [14], which 
emphasizes close relationship to classical catastrophe theory [15]. 

This paper is divided as follows: in section 2, we introduce the general case of the examined 
CPT-symmetric SSH model. In section 3, we link its special case cos# = 0 to a discretized 
Robin square well, and construct a complete set of pseudometrics. In section 4, we repeat this 
procedure for a dual SSH nearest-neighbor interaction. In section 5, we consider the general case 
H = Hssh + 77/, and address the questions of pseudo metric construction and cutoff emergence. 
Section 6 is devoted to discussion and remarks. 

2 The TT-symmetric SSH model 

Inspired by [2], we examine a SSH model with open boundary conditions, complemented by 
the IPT-symmetric boundary term as in eq. 2. Throughout this paper, we set with no loss of 
generality t = 1, A = 1, and denote A = cos# in eq. 1. The coupling constant 7 was taken to 
be purely imaginary in [2], while we assume a more general case 7 = p + \u e C. In order to 
apply powerful tools of linear algebra, we work in a matrix representation of the creation and 
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annihilation operators. The basis of a corresponding two-dimensional Hilbert space is chosen, 
such that 


a 10) = 0, a f |0) = |1) 

a|l) = | 0 ), a + |l) = 0 


In such case, the model from eq. 1 becomes a family of n x n matrices with n = 27V. We 
adopt a strategy of printing the matrices explicitly for the n = 4, as long as the extrapolation 
pattern for higher n G N is clear enough. That means we may write 


h-(4) - hH) + rd 4 ) 

n ~ n SSH ^ n l 


7 —1 — A 

—1 — A —1 + A 

—1 + A —1 — A 

-1 - A 7 * 


(7) 


In parallel with numerical experiments of [2], we illustrate the effects of general complex 
coupling in fig. I. Clearly, there is no reason to confine attention to purely imaginary 7 in our 
IPT-symmetric considerations, as for p reasonably small, there always exists a nonempty interval 
of u) with real spectrum. Size of such interval shrinks with growing p, and finally vanishes 
completely, in the present case for p & 1. 



Figure I: Re [a(H)\ and Im [a(H)\ for the n = 50 SSH model as a function of 9. The values of 
7 are 0.57, 0.57 + 1 and 0.57 + 2. 



Figure II: Exceptional points of eq. 8 for n = 3,4,5 in the (cc,/3) and (p,u}) coordinates. 
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3 The Robin square well 

The construction of 0 for infinite-dimensional quasi-hermitian models is in general a highly 
nontrivial task, which has to be approached perturbatively. An exceptional operator admitting 
exact metric construction (Laplacian on a real interval with imaginary Robin boundary condi¬ 
tions) was introduced in [16]. This inspired a subsequent paper [17], which applied discretization 
techniques on this continuous model, resulting in a TT-symmetric family of matrices 

with 7 = --—— ( 8 ) 

1 — a — ip 

0. Throughout this paper, we shall use the natural 
The coordinate transformation connecting these 


rd 4 ) _ 
^ DR ~ 


7 -1 
-1 -1 
-1 -1 

-1 7* 


which clearly coincides with eq. 7 for A = 
parametrization 7 = p + iuj instead of eq. 8 . 
parametrizations is 


U (1 — a) 2 + f3 2 ’ P (1-a) 2 + /3 2 U 

The domains of quasi-hermicity in both coordinate systems are shown in fig. II. The 
non-shrinking behavior of the domains in the limit n —> 00 agrees with the existence of a 
quasi-hermitian operator in the continuous limit. Now, recall that any n x n metric may be 
expressed as 0 = Ylk=i K n\ n ) ( n \ [18], where |n) are the eigenvectors of iT. Thus, a general 
metric depends on n arbitrary parameters. Equivalently, we can construct n linearly independent 
hermitian (but not necessarily positive) matrices for a given H, such that 

n 

® = J2 £kVik) ( 10 ) 

k =1 

with £k being constrained by the condition 0 > 0. Such a set of V^ is called a complete set 
of pseudometrics. Although the decomposition from eq. 10 may be in general quite arbitrary, 
our aim is to find the pseudometrics in a form, which admits extrapolation to general n G N. 
Already in [17], the authors appreciated the existence of a particular metric 0 with elements 


II j = k 

Qjk = s —i u(p — icuy-h- 1 j > k ( 11 ) 

[ k u(p + icu)- 7_fc_1 j < k 

which served also as a starting point of our considerations. For sufficiently low dimensions, 
we can construct pseudometrics explicitly with the help of symbolic manipulation on any 
computer algebra system. For the present two-site model, after denoting £ = p — icu, we get 
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■p(X) 


1 —\u - 

-iwf 


i uj 1 

—ic o 

—icu£ 

icu£* i uj 

1 

—ic jJ 

i uS,* 2 iu>£* 

\u 

1 



1 — i uj 

-p(3) — 

1 

1 P 

P 1 

1 


i id 1 



■pi 2 ) 


•p(4) 


1 —i u —icu£ 
1 p 1 —\u 

iw 1 p 1 

ic<;£* icu 1 

1 ' 

1 

1 

1 


( 12 ) 


Note that V ^ coincides with eq. 11, while V ^ realizes a discrete operator of parity, and the 
condition = V^H demonstrates the CPT-symmetry of the present model. The two-site 

case suggests an extrapolation pattern for higher n G N with the fc-th pseudometric having 
only 2(n — k) + 1 nonzero antidiagonals. The nonzero elements of such metrics are given by the 
following table 


-i 


> k 

i w £*(i-i-*0 

j ~i 

Al 

P 


i ~ j 

< k, i + j — k even 

1 

V 


i ~ J 

< k, i + j — k odd 


(13) 


A rigorous proof of the above proposition by double induction n —>• n + 1 and k —>• k + 1 
would be a straightforward, although lengthy, generalization of proposition 2 in [17] (which is 
essentially the proof of n —> n + 1 for k — 1), and therefore may be omitted. The positivity of 
the resulting metric as a function of c, is a generally nontrivial problem, which lies outside of 
the scope of the present discussion. We can however make use of the fact, that the pseudo metric 
'P 1 ' 1 ' 1 is a genuine metric, and treat other £{ using tools of perturbation theory. 


4 The dual-SSH model 


A large subclass of tridiagonal quasi-hermitian matrices has a close connection to the theory of 
orthogonal polynomials [19, 20]. Inspired by these models, we define the dual-SSH (or dSSH) 
model, which, in parallel with the operator-matrix correspondence of eq. 7 has its two-site form 


7 —1 — A 

tt( 4 ) _ — 1 + A —1 — A 

*dssH - -i + A -1 - A 

—1 + A 7 * 


(14) 


Schematic drawings of the both dSSH and the original SSH interaction are shown for N = 4 
in fig. Ill, with black circle denoting the boundary term from eq. 2, and solid, respectively 
dashed lines denoting the —1 — A, respectively — 1 + A interaction 

In this section, we study the special case of eq. 14 with 7 = 0. This matrix family 
indeed belongs to a class of tridiagonal models admitting elegant treatment using theory 
of orthogonal polynomials. Recall that any orthogonal polynomial sequence p n (x ) obeys a 
three-term recurrence relation 
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SSH: 


dSSH: 



Figure III: Schemes of the four-site SSH and dSSH models. 


•EPn^pj') ®n(n+l)Pn+l(*^) T ®nnPn(*^) “1“ ®n(n-l)Pn-l(^) (15) 

with a n ( n +i)fl(n+i)n > 0 and a nn £ M. The coefficients may be understood as elements of 
truncated three-diagonal n x n matrices with characteristic polynomial p n (x). Consequently, 
the eigenvalues of such matrices (the roots of p n ) are real and distinct, which is a sufficient 
condition for matrix quasi-hermicity. Moreover [21], a complete set of pseudometrics for such 
models has a very feasible form, with V^ k ’ having only ( 2 k + 1 ) nonzero diagonals, and may be 
constructed by solving the recurrences 

k— 1 k—1 

'^ <1 {k+j){k) P {k+j){k+ 1 ) = ^j a {k+j)[k+l)^(k)(k+j) ( 16 ) 

j =0 j =0 

The matrix elements of eq. 14 obey the condition a n ( n+ i)a( n +i) n > 0 for any 9 kir. At the 
exceptional points 6 = kn the operator is not diagonalizable despite the reality of its spectrum, 
and its quantum-mechanical interpretation is lost. For 6 kn, we may construct pseudometrics 
using eq. 16, or again with the help of computer-based symbolic manipulations. In either case, 
we arrive at the two-site formulae 


•pi 1 ) 


"+ 


1 


+ 


1 ' 

— 

p( 2 ) — 

1 1 

p(3) — 

— — 

pH) = 

1 

+ 

1 1 

+ + 

1 

— 


1 


— 


1 


(17) 


where we have denoted ± = (1 ± A). These formulae suggest an extrapolation pattern with 
k-th pseudometric having 2k + 1 nonzero diagonals and 2{n — k) + 1 nonzero antidiagonals. 
Moreover, the nonzero elements are arranged in a chessboard-like pattern, in accordance with 
[20]. We may write formula for zero elements as 


■T)(k) 

' ij 



for |i + j — k | odd 

for \i — j | > k or \i + j — n — 11 > n — k 


(18) 


Nonzero elements of are given by the following table. In the left column, we have listed 
for comparison also nonzero pseudometric elements for the classical SSH interaction from eq. 7 
with 7 = 0 . The hermicity of such model is apparent from 0 = / being among the admissible 
metrics. 
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SSH dSSH 

+ 

k even, i odd 

+ 

k odd, i odd 

— 

k even, i even 

— 

k odd, i even 

1 

V 

k odd 

1 

k even 

J 


( 19 ) 


5 The inner products 

After examining two special cases 7 = 0 and A = 0, we are ready to address the model of eq. 7 
in full generality. In order to print the resulting matrices explicitly, we define ± = (1 ± A )w, 
where w appearing in the (ij)-th element of the fc-th pseudometric denotes the (ij) -th element 
of the k- th pseudometric of eq. 12. Explicit construction of metric operators remains feasible 
for sufficiently low dimensions. For the t.wo-site model, we get the following results (we shall 
explain shortly why V x ' ! is missing from the list) 


+ 

1 

+ 

_ 1 


+ r 


1 

+- - - 2 + 

■p(3) _ 

• SSH ~ 

- 1 + 

p( 4 ) _ 

' SSH ~ 

1 

+ - 2 - +- 
1 + +- 

+ 1 - 
1 + 

1 

1 


V, 


( 2 ) 


dSSH 


- 2 - 1 ' 

- 2 - +- + 
- +- + + 2 
1 + + 2 


V, 


( 3 ) 

dSSH 


- 1 
- 1 + 
- 1 + 

1 + 


V. 


( 4 ) 


dSSH 


1 

1 

1 


1 


( 20 ) 


Extrapolation formulae for general nSN are clear from these matrices, again with the k- th 
pseudometric having 2(n — k) +1 nonzero antidiagonals. However, as long as we wish to impose 
the same ansatz for pseudometrics V^ k \ k < n — 2, we discover the emergence of a cutoff value, 
for which the pattern ceases to be valid. This phenomenon can be seen already in the case of 
hepta-antidiagonal pseudometrics. We illustrate this behavior on the following three-site models 




+ 2 - 

+ - 

+ 

1 



_ 3 

_ 2 

— 

1 ' 


+ - 2 

+ - 

+ 2 - 

_ 2 

+ 


_3 

+- 

+ 2 - 

+ - 

+ 

+ 2 - 
+ - 

+ - 
+ 2 - 

+ - 2 
Ai 

A, 

+- 2 

+ 2 - 

+- 

+- 

+ 2 - 

p(3) _ 

1 dSSH — 

1 1 

to CO 

!“! 

+- 2 

A 2 

A 2 

+- 2 

+ 2 - 
+ - 

+ 2 

+ 3 

+ 

_2 

+ 2 - 

+- 

+- 2 



- +- 

+ 2 - 

+- 

+ 3 


1 

+ 

+ - 

+ 2 - 




_ 1 + 

+ 2 

+ 3 




( 21 ) 


where H— 2 stands for +(—) 2 . If we insert these matrices into the equation = V^H, 

we obtain complex solutions for A,, which violate the hermicity condition. Thus, in the general 
case of eq. 7, we are able to build just three pseudometric families and -pO) before 

cutoff emergence. Still, the pattern of matrix elements not approaching infinity as n —> 00 is 
strongly conjectured to be preserved, suggesting that the resulting metrics 0 remain bounded 
in appropriate continuous limit. Such a continuous analogue is well-known for A = 0 [16], but 
the examination of the continuous limit for general 6^0 remains open. 
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6 Discussion 


We have successfully constructed a complete set of pseudometrics for the discrete Robin square 
well, as well as for the hermitian SSH model and its dual counterpart. In the general case 
of eq. 1, we encountered a cutoff preventing the construction of pseudometrics beyond three 
particular families. While this may sufficient in many scenarios, we emphasize two reasons for 
the importance of constructing a complete set of pseudometric. As long as the Hamiltonian is 
not the only dynamical observable, the presence of additional observables A,; imposes additional 
constraints on the metric, in the form of equations 

A»e = eA, (22) 

Furthermore, the physical Hilbert spaces 77generated by different metrics are in general 
not unitarily equivalent. In particular, for parametric models, some metrics may have larger 
domains of positivity than others [22, 23]. The choice of a correct physical metric is then dictated 
by other physical principles of, for example, locality. 





Figure IV: Eigenvalues of the n = 50 SSH, dSSH and the model —1 — A —¥ — 1 + A for 
7 = 0.5i + 0.8. 


Spectral behavior of the discussed models is briefly outlined for a single parameter value 
in fig. IV, where we have, for comparison, shown also the operator obtained from eq. 7 by 
transforming all elements —1 — A into — 1 + A. All three models exhibit nontrivial domains of 
observability occurring for non-zero 7. For non-hermitian operators, it may show rewarding, in 
addition to studying the spectrum, to examine their £-pseudospectrum [24, 25] 


a £ (H) = {A eC | ||(i7 - A) _1 || > e” 1 } (23) 

While the pseudospectrum of a general operator may behave quite arbitrarily, pseudospectra 
of quasi-hermitian operators admit a simple characterization based on the spectral theorem. 
Using the decomposition 0 = fUfl, we may write 


1 

p(A, cr(H)) 


< 11(77 - A) -1 1| < 


p(A,cx(T7)) 


(24) 


therefore the c-pseudospectrum is contained in the ^-neighborhood of the spectrum, with 
k = ||fl|| ||f2 -1 ||. Clearly, the pseudospectrum contains information unavailable from the 
spectrum itself, with the resolvent norm being pronounced near the eigenvalues, which are likely 
to complexify by a small perturbation of parameters. 
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Figure V: Pseudospectra of eq. 8 for n — 8 and a — 1, /3 — 1, respectively a = 1, = 0.7 
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